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Abstract 

A  methodology  is  presented  for  constructing  an  exact  stress  solution  for  a 
multiorthotropic-layered  cylinder  for  which  all  layers  have  equivalent  axial 
strain.  Combined  loadings  of  internal  pressure,  external  pressure,  and  axial 
force  are  included  in  the  solution's  formulation.  The  multilayered  solution  is 
developed  by  applying  the  proper  boundary  conditions  to  a  three-dimensional 
elasticity  solution  for  a  monolayered  orthotropic  cylinder.  Theoretical  results 
for  multilayered  cylinders  with  selected  material  layups  and  loading  conditions 
are  compared  with  finite  element  results  to  verify  the  accuracy  of  the  solu¬ 
tion's  formulation.  The  software  imolementation  of  this  solution  provides  the 
design  engineer  with  a  powerful  tool  for  efficiently  using  composite  materials 
in  cylindrical  pressure  vessels. 


1.  Introduction 

Composite  materials  have  become  important  in  structural  applications  requiring 
high  stiffness  and  low  weight.  This  situation  also  exists  in  a  current  problem 
with  Army  cannon.  The  Army  is  interested  in  using  comDOsite  materials  to 
lengthen  cannon  while  still  maintaining  the  inertial  characteristics  of  the 
shorter  cannon.  The  design  and  analysis  of  these  new  cannon  require  a  stress 
solution  for  thick-walled  composite  cylinders. 


1.1  Monolayered  Stress  Solutions 

The  simplest  type  of  cannon  construction  involves  the  use  of  only  one  material 
layer.  The  well-known  stress  solution  by  Lam6  in  1852,  as  outlined  in  [1],  >s 
sufficient  to  characterize  the  state  of  stress  for  a  monolayered  isotropic 
cylinder,  lekhnitskii  [2]  developed  a  more  general  stress  solution  for  -i  mono¬ 
layered  anisotropic  cylinder.  Recently,  O’Hara  [3]  investigated  lekhn' ;ski i ' s 
equations,  simplifying  them  for  the  case  of  an  orthotropic  material  under 
loadings  of  internal  pressure,  external  pressure,  and  axial  force.  Shaffer  [4] 
and  Bieniek,  Spillers,  and  Freudenthal  [5]  also  investigated  stress  solutions 
for  monolayered  cylinders  made  from  special  types  of  orthotroDic  materials. 


1.2  Multilayered  Stress  Solutions 

In  real  applications,  however,  cylinders  are  often  constructed  from  many  com¬ 
posite  layers  where  the  fibers  are  wound  or  layed  up  at  various  angles.  Most 
commonly,  fibers  are  used  in  +  and  -  wrap  angle  pairs,  where  each  pair  can  be 
viewed  as  a  single  orthotropic  layer.  The  whole  structure  can  be  considered  a 
multiorthotropic-layered  cylinder.  Stress  solutions  for  bi layered  cylinders 
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have  been  investigated  by  Shaffer  [6]  and  Witherell  and  ScavuTIo  [7].  For 
multiorthotropic-layered  cylinders,  stress  solutions  have  been  developed  by 
Bieniek,  Spillers,  and  Freudenthal  [5],  Roy  and  Tsai  [8],  and  Witherell  [9]. 

The  solution  by  Roy  and  Tsai  is  the  most  general  and  is  a  very  good  approximate 
stress  solution  for  a  multiorthotroDic-layered  cylinder  with  generalized  plane- 
strain  boundary  conditions  under  loadings  of  internal  and  external  pressure  and 
axial  force. 


1.3  Exact  Multilayered  Stress  Solution  For  Generalized  Plane-Strain  Boundary 
Conditions 

This  paper  is  similiar  to  the  work  of  Roy  and  Tsai  [8],  but  whereas  the  develoo- 
ment  of  their  multilayered  solution  does  not  include  the  total  contribution  to 
the  hoop  strain  eauivalence  condition  from  the  radial  and  hoop  stress  produced 
by  an  axial  force  component,  the  solution  discussed  in  this  oaoer  does.  The 
inclusion  of  these  two  hooo  strain  contributions  in  the  solution  is  not  trivial. 
Put  it  does  result  in  a  multilayered  stress  solution  that  is  more  rigorously 
correct  than  that  of  Roy  and  Tsai.  For  several  material  layups  and  loading  con¬ 
ditions  the  author  has  investigated,  the  difference  between  the  stress  solution 
in  this  oaoer  and  that  of  Roy  and  Tsai  is  small  (less  than  10  percent).  To  some 
degree  this  small  difference  is  to  be  expected,  since  the  components  neglected 
in  the  Roy-Tsai  solution  are  of  a  secondary  nature  being  a  conseauence  of  dis¬ 
similar  Poisson  contraction  or  expansion  among  axially  loaded  layers.  However, 
there  may  oe  some  material  layuos  and  loading  conditions  where  this  difference 
is  greater.  In  addition,  the  exact  solution  presented  in  this  Daper  provides  a 
solid  foundation  for  developing  other  stress  solutions  and  eliminates  any  unnec¬ 
essary  uncertainties  that  would  be  associated  with  building  on  a  less  accurate 
approximate  stress  solution.  For  example,  the  solution  that  follows  is  a  neces¬ 
sary  part  of  the  complete  formulation  for  the  exact  thermal  stress  solution  for 
a  multilayered  cylinder  with  generalized  plane-strain  boundary  conditions.  This 
oaoer  describes  the  methodology  of  constructing  the  exact  elastic  stress  solu¬ 
tion  for  a  multiorthotropic-layered  cylinder  under  loadings  of  internal  ana 
external  pressures  and  axial  force.  A  comparison  is  then  made  between  the  exact 
multilayered  solution  and  a  finite  element  solution  to  verify  the  solution  pro¬ 
cedure. 


2.  Geometry,  Material,  and  Loading  Definition 

A  multilayered  cylinder  can  be  viewed  as  an  assembly  of  many  single- layered 
cylinders.  It  is  fitting,  therefore,  to  begin  with  a  review  of  the  monolayered 
orthotropic  cylinder  problem.  There  are  actually  two  monolayered  solutions  that 
are  used  to  construct  the  generalized  plane-strain  monolayered  stress  solution 
for  internal  and  external  pressures  and  axial  force.  The  first  solution,  called 
SOLI,  is  for  internal  and  external  pressure  loadings  with  plane-strain  boundary 
conditions  (zero  axial  strain).  The  second  solution,  called  S0L2,  is  for  axial 
loading  with  generalized  plane-strain  boundary  conditions  (constant  axial 
strain) . 


2.1  Mono layered  Problem 

In  the  monolayered  case,  the  cylinder  is  assumed  to  be  long  with  ends  that  are 
either  fixed,  as  in  SOLI,  or  parallel  such  that  the  axial  strain  is  uniform,  as 
in  S0L2 .  The  geometry  and  loading  for  SOLI  and  S0L2  are  shown  in  Fig.  l.  The 
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cylinder  has  an  inner  radius  'a'  and  an  outer  radius  'b'.  For  SOLI,  the 
cylinder  is  subjected  to  internal  pressure  ' p '  and  external  pressure  'q',  and  it 
has  an  axial  force  Fzj  necessary  to  enforce  the  plane-strain  axial  constraint. 
For  S0L2,  the  cylinder  is  only  subjected  to  an  axial  force  fZ2-  For  bo*1"1  SOLI 
and  S0L2,  a  cyl indrical ly-orthotropic  material  is  assumed  with  its  principal 
axes  coincident  with  the  cylindrical  coordinate  system  defining  the  cylinder 
geometry.  An  orthotropic  material  is  characterized  by  nine  independent  material 
constants  consisting  of  three  engineering  moduli i  (E1.E2.E3),  three  shear  modu- 
lii  ( S12 > S23 » S31 ) ,  and  three  Poisson's  ratios  ( * ^23 * y31 ) •  The  numbers  1,2,3 
indicate  the  principal  material  directions.  For  the  above  assumptions,  the 
1,2,3  directions  correspond  to  the  radial,  hoop,  and  axial  directions  of  the 
cylinder  (r,0,z).  In  addition,  since  the  principal  stress  directions  correspond 
to  the  principal  directions  of  both  the  cylinder  geometry  and  the  applied 
loadings,  shear  effects  are  eliminated,  and  the  number  of  material  constants 
necessary  for  the  analysis  reduces  to  six. 


3.  Stress  Equations  For  a  Monolayered  Orthotropic  Cylinder 

..ekhnitski  i '  s  two  solutions  for  a  cylinder  with  one  anisotroDic  layer  were 
simplified  by  O'Hara  [3]  for  the  case  of  a  cyl indrical ly-orthotrooic  material. 
The  equations  for  these  two  solutions  (SOLI  and  S0L2)  are  given  in  the  Appendix 
and  are  identical  to  those  found  in  O'Hara's  report  with  several  corrections  for 
some  minor  typographical  errors.  SOLI  and  S0L2  each  include  three  stress 
eauations  corresponding  to  the  r,  0,  and  z  directions.  SOLI  also  contains  an 
equation  which  defines  Fz^,  the  axial  force  necessary  to  enforce  the  plane- 
strain  axial  constraint.  S0L2  also  has  an  additional  eauation  that  contains  a 
complicated  constant  (T)  used  in  the  three  other  stress  equations. 


3.1  Monolayered  Stress  Components  at  Inner  and  Outer  Surfaces 

In  developing  the  multilayered  solution,  applying  the  correct  boundary  con¬ 
ditions  at  the  interface  of  two  orthotropic  layers  necessitates  the  use  of 
stress  values  evaluated  at  the  inner  (r=a)  and  outer  (r=b)  surfaces  of  each 
layer.  This  evaluation  process  leads  to  six  stress  equations  for  SOLI  and  six 
stress  equations  for  S0L2  given  below. 

SOLI  -  Stress  equations  evaluated  at  inner  radius  (r  =  a): 


ffr,a  =  RAP»o  +  RAQ-o 
RAP  =  -1  ,  RAQ  =  0 

a0,a  =  TAP*P  +  TAQ*q 

k(l+C02k)  -2kC0k"1 

TAP  = - |r-  .  TAQ  =  -- ?--- 

(1-Co  )  (1-Co  1 

az  a  =  ZAP* 0  ♦  ZAQ*q 


(Ai  3-A23  •  TAP )  -A'j-i^TAQ 

ZAP  =  — .  ,  ZAQ  =  — -t - 


*33 


*33 


(2) 
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SOLI  -  Stress  equations  evaluated  at  outer  radius  (r  »  b): 

crr  ^  =  RBP*p  +  RBQ*q  (4) 

RBP  =  0  ,  RBQ  =  -1 

a9,b  =  TBP*p  +  TBQ.q  (5) 

TBP  =  -TAQ*C02  ,  TBQ  =  -TAP 

az,b  s  ZBP*p  +  ZBQ»q  (6) 

ZBP  =  -ZAQ*C02 

The  axial  force  equation  for  SOLI  can  also  be  rewritten  in  a  similiar  form 

Fzl  =  FZP-P  +  FZQ*q  (7) 


2Ai  3 

ZBQ  =  -ZAP  +  ---- 
a33 


FZP  = 


2naf 

[ (TAP+TAQ+1 ) (A23-A13)] 

a33  ' 

(1-k2) 

-  A 


FZQ  = 


2?rb2 

[ (TBP+TBQ+1 ) ( A23~A13 ) ] 

a33 

(1-k2) 

23 


'23 


wnere 


C°  =  b 


and  the  comDonents  A.jj  are  the  elements  of  the  compliance  matrix  for  an 
orthotrooic  material,  as  given  in  the  generalized  Hooke's  Law, 


[€]  =  [A]  [o] 


' k ’  is  an  orthotropic  material  constant  given  by 


k  ’  511 

_  ^22 

where,  in  general, 


Si  i  =  A 


■>  J 


Ai3_A3j 

A33 


S0L2  -  Stress  equations  evaluated  at  inner  radius  (r  =  a): 


°r , a  =  0 
CT9.a  =  TAF.Fz2 
TAF  =  !p  ( 1+TAP+TAQ ) 

°z,a  ~  ZAF*Fz2 


ZAF  = 


1 

T 


A23 • TAF 
a33 


(8) 


(9) 


(10) 


(11) 

(12) 


(13) 
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(14) 


S0L2  -  Stress  equations  evaluated  at  outer  radius  (r  =  b): 


ar , b  =  0 
a9 , b  =  TBF«Fz2 
TBF  =  ^  ( 1+TBP+TBQ) 
a2,b  =  ZBF*Fz2 


where 


and 


ZBF 


l  A23  »TBF 

^  a33 


fill  -  fl22 


T  =  it(  b2-a2 ) 


27rh  b2-a2 

A33  2 


(a13+a23> 


b2  ( 1-C0k+1 ) 2 

d-c02k) 


a2(l-C0k  1)2  (Ai3-kA23) 

d-c02k)  (l-k) 


(Ai3+kA23 ) 

(l+k) 


(15) 


(16) 


(IT) 


( 18) 


In  each  of  the  twelve  stress  equations  above,  there  are  constants  (two  for  SOLI 
and  one  for  S0L2),  e.g.,  TAP  and  TAQ  for  the  ag<a  stress  of  SOLI  and  TBF  for  the 
oq t b  of  S0L2 ,  which  determine  the  magnitude  of  the  contribution  to  the  given 
stress  caused  by  either  the  internal  or  external  pressure  for  SOLI  or  the  axial 
*or ce  for  SOL2.  These  constants  are  material-  and  geometry-dependent,  ana  eacn 
nas  a  three-letter  name.  The  first  letter  corresponds  to  the  stress  direction 
(r,9,z),  the  second  letter  signifies  the  point  at  which  the  stress  has  been 
evaluated  (a,b),  and  the  third  letter  corresponds  to  either  the  pressure  (q,p) 
it  is  applied  to  or  the  axial  force  Fz2.  For  example,  TAQ  is  the  constant  for 
the  Theta  stress,  evaluated  at  r  =  a  (A),  magnifying  the  external  pressure  a 
(Q),  and  ZBF  is  the  constant  for  the  Z  stress,  evaluated  at  r=b  ( B ) ,  magnifying 
the  axial  force  Fz2. 


4.  Stress  Solution  For  a  Multiorthotropic-Layered  Cylinder 

In  this  section  we  discuss  the  methodology  of  constructing  the  stress  solution 
for  a  multilayered  cylinder,  where  each  layer  has  material  and  geometry  defini¬ 
tions  identical  to  those  previously  discussed  for  the  monolayered  case.  The 
solutions  SOLI  and  S0L2  are  used  as  the  basis  for  the  multilayered  solution. 


4.1  Decomposition  of  the  Generalized  Plane-Strain  Multilayered  Solution 

The  complete  multilayered  solution  for  internal  pressure,  external  pressure,  and 
axial  force  loadings  with  generalized  plane-strain  boundary  conditions  is 
constructed  from  the  superposition  of  three  separate  multilayered  solutions. 
Figure  2  shows  the  decomposition  of  the  overall  problem  into  its  three  parts. 

The  first  solution,  called  MS0L1,  is  the  multilayered  plane-strain  solution  for 
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internal  and  external  pressure.  This  solution  requires  finding  the  set  of 
interface  pressures  such  that  the  hoop  strain  of  each  layer  at  every  two- layer 
interface  is  equal.  Stated  another  way,  the  radial  displacement  must  be  con¬ 
tinuous  through  the  cylinder  (no  gaps).  This  solution,  which  uses  SOLI,  has 
already  been  obtained  by  Witherell  [9].  The  total  axial  force  necessary  to 
enforce  the  plane-strain  boundary  conditions  for  MS0L1  is  denoted  as  Fzit  and  is 
the  summation  of  the  axial  forces  for  each  of  the  layers  (Fz^)  as  given*in  El, . 
(7).  The  second  solution,  called  MS0L2  (or  the  axial  force  solution),  is  the 
multilayered  solution  obtained  by  applying  an  axial  force  to  each  layer  such 
that  all  the  layers  have  equivalent  axial  strain.  For  this  solution,  each  layer 
is  considered  independent,  and  thus  gaps  or  overlaps  result  at  each  layer  inter¬ 
face  once  the  axial  forces  are  applied.  These  gaps  or  overlaps  are  taken  care 
of  in  the  third  solution.  It  should  be  noted  that  the  axial  strain  produced  in 
each  of  these  layers  is  the  total  axial  strain  for  the  overall  solution.  The 
total  axial  force  that  results  for  the  MS0L2  solution  is  equal  to  the  sum  of  the 
axial  forces  on  each  layer  (FZ2>  and  is  denoted  as  ?z2t4  Finally,  the  third 
solution,  called  MS0L3,  is  a  plane-strain  solution  for  a  set  of  interface 
pressures  that  eliminates  the  gaps  or  overlaps  produced  by  MS0L2.  MS0L3  is 
called  the  plane-strain  equilibrating  solution  because  it  equiliorates  the  hoop 
strain  at  the  layer  interfaces  by  removing  the  gaps  or  overlaps  produced  by 
MS0L2.  The  form  of  the  solution  for  MS0L3  is  identical  to  MS0L1,  except  that  it 
does  not  have  internal  or  external  pressures  acting  on  the  overall  cylinder 
(this  was  accomplished  in  MS0L1),  only  pressures  acting  on  the  interfaces  of 
each  layer.  Also,  since  MS0L3  is  a  plane-strain  solution,  it  does  not  affect 
the  total  axial  strain  of  the  MS0L2  solution.  The  MS0L3  solution  does,  however, 
produce  a  total  axial  force  which  must  be  accounted  for.  The  total  axial  force 
for  MS0L3  is  denoted  as  FZ3t  and  is  obtained  by  the  summation  of  all  the  layer 
axial  forces  ( F Z3 )  (FZ3  is  of  the  same  form  as  Fz^  in  Eq.(7))  produced  by  the 
interface  pressures  of  MS0L3.  The  net  axial  force  for  the  overall  solution, 
prescribed  for  a  given  problem,  is  denoted  as  Fnet  and  "is  eciual  to  the  sum  of 
P2lt'  Fz2t'  and  pz3t*  The  following  section  contains  the  necessary  aetails  for 
constructing  the  MS0L1,  MS0L2,  and  MS0L3,  as  well  as  combining  these  three  solu¬ 
tions  to  obtain  the  overall  solution. 


4.2  The  MS0L1  Multilayered  Solution 

MS0L1  is  the  solution  for  the  multilayered  cylinder  with  internal  and  external 
pressure  under  plane-strain  boundary  conditions  and,  as  mentioned  earlier,  has 
already  been  developed.  The  following  discussion  gives  the  important  aspects  of 
this  solution  also  used  to  develop  MS0L3 . 


4.2.1  Step  1  For  MS0L1 

the  first  step  in  constructing  MS0L1  is  to  equate  the  ci rcumf erenti al  (hoop) 
strain  of  adjacent  orthotropic  layers  at  their  interface.  In  other  words,  we 
eauate  the  hoop  strain  at  r  =  b  of  the  ith  layer  with  the  hooo  strain  at  r  =  a 
of  the  i+l  layer 


se.bt1)  =  e9  ,  a ( i+l )  (19) 

Recalling  that  [e]  =  [A]  [o] ,  1,2,3  correspond  to  the  r,e,z  directions,  and  by 
using  the  stress  equations  of  SOLI  with  q(i-l)  =  p(i),  we  have 

G  i  1  *  ( i  - 1 )  *  G-j2*q(i)  +  Gi3*q(i+1)  =  0  (20) 
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where 


Gii  =  "Co2  ( i )  *622^  )  •  TAQ ( i ) 

G-j  2  =  “  [012  ( 1  )  "  J3i2(i  +  1)  +  022  ( 1 )  ‘TAP(  i  )  +  022  ( i  +  l )  *TAP  ( i  +  1 )  ] 

Gi'3  =  “022^  +  l )  *TAQ(i+l) 


4.2.2  System  of  Equations  For  MS0L1 

For  each  two-layer  combination  there  is  one  equation  defining  the  hoop  strain 
equivalence  condition.  For  a  cylinder  with  'N'  orthotropic  layers,  there  are 
N-l  equations  and  N-l  unknowns.  Setting  up  these  equations  in  matrix  form  and 
noting  that  p(l)  and  q(N)  are  the  prescribed  internal  and  external  pressures  of 
the  overall  cylinder  results  in 


r  1 

— 

—  — 

I 

G12  g13 

l 

0(1) 

*Gn*p(  1 ) 

|  G21  g22  g23  0 

q  ( 2 ) 

0 

j  g31  g32  g33 

i  _  _ 

q  ( 3 ) 

0 

0 

- 

’ 

|  gN-2 , 1  gN-2,2  gN-2 , 3 

1  gN-1,1  gN-1 , 2 

q ( N-l ) 

-Gn_i f  3  *d(N ) 

1_  *  — 

_  — 

L  J 

or 

[G]  [Q]  =  [R]  (22) 


4.2.3  Stress  Distribution  in  Each  Layer  For  MS0L1 

This  system  of  equations  can  be  easily  solved  to  determine  the  interface 
pressure  vector  [Q] .  By  recalling  that  q(i-l)  =  p(i),  we  can  now  calculate  the 
stress  distribution  in  layer  i  using  the  general  stress  equations  for  SOLI  found 
in  the  Appendix  with  p(i)  and  q(i)  as  input. 


4.2.4  Total  Axial  Force  For  MS0L1 

In  addition  to  finding  the  stress  distribution  in  ^yer  i,  the  axial  force 
necessary  to  constrain  layer  i  to  zero  axial  strain  can  be  calculated  using 
Eq.  (7)  with  p ( i )  and  q ( i )  as  input.  By  summing  the  axial  forces  on  each  layer, 
we  can  determine  the  total  axial  force,  Fzjt,  for  MS0L1 

N 

Fzn  »  I  Fzim 

i  =  1 


(23) 


4.3  The  MS0L2  and  MS0L3  Multilayered  Solutions 

Next,  we  proceed  with  the  MS0L2  and  MS0L3  solutions,  which  need  to  be  looked  at 
together.  Since  the  internal  and  external  pressure  loadings  were  already  taken 
care  of  in  MS0L1,  they  are  not  included  here.  However,  we  need  to  bring  the 
the  total  axial  force  F^lt  MS0L1  into  this  discussion  because  it  is  part  of 
the  net  axial  force  on  the  cylinder. 


4.3.1  Step  1:  Equating  the  Axial  Strain  of  Each  Layer  (Each  Layer  Independent) 

We  accomplish  this  by  eauating  the  axial  strain  at  the  external  surface  of  each 
layer  (r=b)  to  that  of  layer  1 


, b ! 1 )  =  £z,p( i )  (24) 

By  substituting  Hooke's  Law  with  the  stress  components  from  S0L2,  we  arrivi  at 
an  exDression  for  the  axial  force  of  layer  i  in  terms  of  the  axial  force  of 
layer  l  and  the  material  and  geometry  properties  of  layer  l  and  layer  i 


wnere 


_  ...  _  , , ,  CAX ( 1 ) 

Fz2(i)  =  F22<1)  CAXli) 


(  25) 


CAX(i)  =  A32 ( i ) • TBF ( i )  +  A33 ( i ) • ZBF ( i )  (26) 

4.3.2  Step  2:  Equating  the  Hoop  Strain  at  Each  Layer  Interface 

This  steD  reauires  consideration  of  only  MS0L2  and  MS0L3  solutions,  since  the 
nooD  strain  eauivalence  at  each  layer  interface  was  already  satisfied  n  MS0L1 

ee,p(i)  =  ee,a(^+l)  (191 

The  hooD  strain  at  the  interface  of  layer  i  and  i+1  has  contributions  from  the 
two  stress  solutions  MS0L2  and  MS0L3: 

£0 ,b ( i )  +  ee.bi1)  =  £0,a(i+l)  +  ee,a(i+l)  (27) 


Axial  Force 
MS0L2 


Plane-Strain  Equilibrating 
MS0L3 


4.3.2. 1  Axial  Force  Contribution  to  ee.bM)  and  £0,a(i+ll 

Jsmg  Hooke's  Law  and  the  S0L2  solution  gives 

£0,b(i)  =  F  z  2 ( i ) *CTB( i ) 

where 

CTB(i)  =  A22 ( i ) *TBF ( i )  +  A23 ( i ) • ZBF { i ) 

Simi 1 iarly , 


(28) 


(29) 
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e9>a(i+l)  *  Fz2(i+l)*CTA(i+l) 

where 

CTA(  1+1 )  =  A22  ( i+l )  *TAF ( i  +  1 )  +  A23  ( i  +  1 )  *ZAF(  i+.l ) 

4. 3. 2. 2  Plane-Strain  Equilibrating  Contribution  to  eg^(i)  anc^  e9,a('i+l) 
From  an  earlier  derivation  of  the  plane-strain  solution  MS0L1 

e© ,  b  ( -i )  =  ee,a('i  +  l) 

reduces  to 

Gil*P(i“l)  +  Gi2*<l('')  +  Qi3*cl(''+1)  =  0 


(30) 

(31) 


(19) 


(20) 


4. 3. 2. 3  Combining  the  Axial  Force  and  Plane-Strain  Equilibrating  Contributions 
Into  Eq.  (27) 

Fz2(  i  ) -CTB(  i  )  -  Fz2(i+l)»CTA(i+l)  +  G1-1*q(i-l)  «•  Gi2*q(i)  +  G1-3*q(i+l)  =  0  (32) 

Now  substituting  Eq.  (25)  for  Fz2(i)  and  Fz2(i+1)  results  in 


=22(1) -CAX(l) 


CTB  ( i  ) 

CTA( i+1 ) 

CAX ( i  ) 

CAX (i+1) j 

+  G-ji*q(i-l)  +  Gi2*a('i)  +  Gi3*a(i  +  1)  =  0 


At  this  point,  the  only  unknowns  are  the  axial  force  Fz2(l)  and  the  interface 
Dressures  q(i-l),  q(i),  and  q(i+l).  Step  3  enables  us  to  get  Fz2(i)  in  terms  of 
the  interface  pressures  and  other  known  quantities,  and  thus  we  can  solve  the 
entire  system  of  equations. 


4.3.3  Step  3:  Summation  of  Axial  Forces 

This  steo  considers  the  total  axial  force  contributions  from  MS0L1,  MS0L2,  ana 
MSQL3 


Fnet  =  Fzlt  +  Fz2t  +  Fz3t 

Both  Fnet  and  Fzjt  are  known  values  (Fzjt  was  calculated  in  MS0L1 ) 

N  N  . 

Fz2t  =  I  fz2  ( 1  )  =  Fz2(l).CAX(l).(  I  6^7-T) 


i  =1 


i  =  l 


Also 


N  N 

Fz3t  =  l  fz3^)  =  l  (FZP(i).p(i)  +  FZQ  ( i  )  *q  (  i  )  } 
i =1  i=l 

For  the  equilibrating  solution,  p(l)  and  q(N)  =  0  and  p(i+l)  =  q(i 


!  34 ) 


(35) 


(36) 


N-l 

Fz3t  =  l  (FZP(i+l)  +  FZQ( i ) ) *q( i )  (37) 

i  =1 
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N  .  N-l 

Fnet  *  Fzlt  ♦  Fz2(l).CAX(l).(  I  5™Ty)  +  X 


i=l 


i=l 


(FZP( i+1 )  +  FZQ( i ) ) *q ( i ) 


(38) 


Solving  for  FZ2<1)  gives 


N-l 

‘net  “  Fzlt  ~  X 
i=l 


(FZP( i+1 )+FZQ( i ) ) *q( i ) 


Fz2(D  = 


(39! 


CAX ( 1 ) • (  I 


CAX(i) 


Now,  back  substituting  FZ2(1)  into  the  hoop  strain  equivalence  condition  for 
interface  ( i )  gives 


r 


CEQ(  i 


N-l 

'Fnet  +  Fzlt  f  X 

j-1  L 


i 


( FZP ( j+l ) +FZQ ( j ) ) *0 ( j 


+  G1-j*q(i-l)  +  G  t  2  *  d  ( i )  +  G-j3*q(i  +  l)  =  0  (40) 

where  CEQ(i)  is  given  by 

CTB(i)  _  CTAM+1) 

CAxIi)  ‘  CAXIi+lI 

CEQ(i)  - - -  Ml; 

N  . 

CAX(j)J 

4.3.4  Step  4:  Solution  Assembly 

In  matrix  form,  the  set  of  N-l  hoop  strain  equivalence  equations  is  given  by 

([F]  +  [G] )  [QQ]  =  FZN.[C]  (42) 

where  the  components  of  [F]  are  given  by 

Fij  =  CEQ ( i ) • ( FZP ( j+l ) +FZQ( j ) )  .  (1  <  i  ,  j  <  N-l)  (43) 

and  the  components  of  the  matrix  [G]  and  the  vector  [QQ]  are  the  same,  as  given 
by  [G]  and  [Q]  in  Eqs.  (21)  and  (22),  The  constant  FZN  is  given  by 

FZN  =  Fnet  -  Fzlt  (44) 

and  the  components  of  the  vector  [C]  are  given  by 

C-j  =  CEQ(i)  ,  (1  *  i  S  N-l)  (45) 

The  above  matrix  equation  can  be  easily  solved  to  determine  the  components  of 
the  vector  [QQ] ,  the  interface  pressures. 
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4.4  Stress  Distribution  in  Each  Layer  For  the  Generalized  Plane-Strain 
Multilayered  Stress  Solution 

In  order  to  obtain  the  stress  distribution  in  layer  i  for  the  generalized  plane- 
strain  multilayered  solution,  the  following  procedure  must  be  employed: 

1.  The  total  interface  pressure  vector  [QT]  must  be  obtained  by  adding  [Q] 
found  from  Eq.  (22)  with  [QQ]  found  from  Eq.  (42) 

CQT]  =  [Q]  +  [QQ]  (46) 

2.  The  axial  force  components  FZ2(i)  must  be  obtained  by 

•  calculating  Fz2(l)  (Eq.  (39))  using  the  components  of  [QQ]  for  q(i)  and 
using  Fz^t  obtained  using  [Q]  from  MS0L1  in  connection  with  Eqs.  (7)  and 

(23); 

•  using  Fz2(i)  and  Eq.  (25),  the  Fz2(i)  components  can  now  be  obtained 
directly. 

The  stress  distribution  in  layer  i  may  now  be  obtained  by  superimposing  the  SOLI 
solution  using  the  components  of  [QT]  as  input  for  q ( i )  and  p(i)  with  the  S0L2 
solution  using  Fz2(i)  as  input  for  Fz2. 

5.  Solution  Verification 

The  exact  multilayered  solution  was  used  to  investigate  the  stress  and  strain 
distribution  within  a  ten-layered  cylinder  (N=10).  The  cylinder  had  an  insioe 
radius  a(l)  =  1  inch,  an  outside  radius  b(10)  =  2  inches,  and  each  of  the  ten 
layers  was  0.1  inch  thick.  The  material  used  was  an  IM6/eooxy  with  a  55  Dercent 
fiber-volume  ratio.  The  layuo  considered  consisted  of  five  hooD-axial  Dairs 
starting  at  the  inside  radius  of  the  cylinder.  The  material  properties  for  the 
hoop  and  axial  fiber  orientations  are  given  in  Table  1.  The  loading  conditions 
for  the  case  considered  included  internal  pressure  (p(l)  =  2  psi ) ,  external 
Dressure  ( q ( 10 )  =  1  psi),  and  axial  force  (Fnet  =  10  lbs).  The  results  for  this 
case  are  shown  in  the  four  plots  of  Fig.  3.  The  plots  show  the  r,e,z  stress  and 
strain  distributions  as  a  function  of  radial  position  in  the  cylinder. 


5.1  Finite  Element  Solution  Used  in  Comparison 

The  multilayered  solution  for  the  above-mentioned  case  was  also  compared  to  a 
finite  element  solution  for  the  same  case.  The  ABAQUS  finite  element  code  was 
used  to  produce  the  stress  results  used  in  the  comparison.  The  finite  element 
model  contained  20  eight-node  quadratic  axisymmetric  elements  with  two  equal¬ 
sized  elements  used  to  model  each  of  the  ten  orthotrooic  layers.  The  comoanson 
was  limited  to  the  radial,  hoop,  and  axial  stress  values  evaluated  at  the  inner 
and  outer  radii  of  each  orthotropic  layer.  Table  2  contains  the  results  of  this 
comparison . 
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6.  Discussion  of  Results 


6.1  Comparison  of  Exact  Solution  With  Finite  Element  Solution 

A  comparison  of  the  theoretical  solution  with  the  finite  element  solution  for 
the  radial,  hooD,  and  axial  stresses  at  the  inner  and  outer  surfaces  of  each 
orthotropic  layer  is  shown  in  Table  1.  For  all  three  stress  components,  the 
comparison  shows  excellent  agreement  between  the  results  generated  by  the  theory 
and  those  from  the  finite  element  solution.  For  the  finite  element  solution, 
the  stresses  were  obtained  at  the  nodes  and  were  somewhat  less  accurate  than  the 
stresses  obtained  at  the  integration  points  of  the  element.  This  error, 
although  very  small,  is  seen  most  clearly  by  the  fact  that  the  radial  stress 
value  at  the  inner  surface  of  the  overall  cylinder  is  not  equal  to  the 
prescribed  value  of  -2  psi. 


6.2  Stress  and  Strain  Distributions  Produced  by  Exact  Solution 

In  Fig.  3,  the  stress  and  strain  distributions  are  displayed  for  the  multi¬ 
layered  cylinder  with  the  same  geometry  and  loading  conditions  as  the  case 
already  discussed.  As  can  be  seen  in  the  blot  containing  the  radial  and  axial 
st-ess,  the  radial  stress  varies  continuously  from  -2  psi  at  the  inside  radius 
to  -1  psi  at  the  outside  radius.  As  is  also  seen  in  the  plots,  the  axial  layers 
take  uo  the  majority  of  the  axial  load,  and  the  hoop  layers  sustain  the  highest 
levels  of  hoop  stress.  If  the  axial  stress  was  integrated  over  the  end  surface 
area,  it  would  be  equal  to  the  net  axial  force  prescribed  for  the  cylinder 
(Fnet-1°  lbs).  The  strain  plots  show  the  expected  smooth  and  continuous  curve 
for  the  hoop  strain  and  the  uniform  axial  strain,  which  the  generalized  plane- 
strain  boundary  conditions  require. 


7.  Conclusion 

An  exact  generalized  plane-strain  elastic  stress  solution  for  a 
multiorthotrobic-layered  cylinder  has  been  developed  for  loadings  of  interna' 
pressure,  external  pressure,  and  axial  force.  The  software  implementation  of 
this  stress  solution  has  been  shown  to  be  in  excellent  agreement  with  finite 
element  results. 
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Table  1.  Material  Properties  For  IM6/Epoxy  55%  Fiber-Volume  Ratio 


Fiber 

Direction 

Er  (Mpsi) 

E e  (Mpsi) 

Ez  (Mpsi) 

ur9 

— 

y0Z 

j 

yzr  | 

Hoop 

1 

Hpsj|||i 

0.0152 

0.3147 

0.3991  ! 

Axial 

wmmm 

0.3991 

0.0152 

0.3147  1 
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Table  2.  Comparison  Between  Multilayered  Solution  and  ABAQUS  Finite  Element  Solution 

Generalized  Plane-Strain  Boundary  Conditions 
Using  20  CAX8  Elements 
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228  2.229 


Fig.  1  Monolayered  cylinder  geometry  and  loading  for  SOLI  and  S0L2. 
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RADIAL  POSITION  (IN.)  RADIAL  POSITION 


Appendix 


SOLI  -  Plane-Strain  Solution  For  a  Monolayered  Orthotropic  Cylinder  Under 
Internal  Pressure  p  and  External  Pressure  q 
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where  'r'  -is  the  radial  position  in  the  cylinder. 


S0L2  -  Generalized  Plane-Strain  Solution  For  a  Monolayered  Orthotropic  Cylinder 
Under  Axial  Force  Fz2 
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TECHNICAL  REPORT  INTERNAL  DISTRIBUTION  LIST 


NO.  OF 
COPIES 


CHIEF,  DEVELOPMENT  ENGINEERING  OIVISION 

ATTN:  SMCAR-CCB-DA  1 

-DC  1 

-01  1 

-DR  1 

-OS  (SYSTEMS)  1 

CHIEF,  ENGINEERING  SUPPORT  DIVISION 

ATTN:  SMCAR-CCB-S  1 

-SO  1 

-SE  1 

CHIEF,  RESEARCH  OIVISION 

ATTN:  SMCAR-CCB-R  2 

-RA  1 

-RE  1 

-RM  1 

-RP  1 

-RT  1 

TECHNICAL  LIBRARY  5 

ATTN:  SMCAR-CCB-TL 

TECHNICAL  PUBLICATIONS  &  EDITING  SECTION  3 

ATTN:  SMCAR-CCB-TL 


OPERATIONS  DIRECTORATE 
ATTN:  SMCWV-ODP-P 

DIRECTOR,  PROCUREMENT  DIRECTORATE  1 

ATTN:  SMCWV-PP 

DIRECTOR,  PRODUCT  ASSURANCE  DIRECTORATE  1 

ATTN :  SMCWV-QA 


NOTE:  PLEASE  NOTIFY  DIRECTOR,  BENET  LABORATORIES,  ATTN:  SMCAR-CCB-TL,  OF 

ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST 


NO.  OF 
COPIES 

ASST  SEC  OF  THE  ARMY 

RESEARCH  AND  DEVELOPMENT 

ATTN:  DEPT  FOR  SCI  AND  TECH  1 

THE  PENTAGON 

WASHINGTON,  D.C.  20310-0103 
ADMINISTRATOR 

DEFENSE  TECHNICAL  INFO  CENTER  12 

ATTN:  OTIC-FDAC 
CAMERON  STATION 
ALEXANDRIA,  VA  22304-6145 


COMMANDER 
US  ARMY  ARDEC 

ATTN:  SMCAR-AEE  1 

SMCAR-AES ,  BLDG.  321  1 

SMCAR-AET-0 ,  BLDG.  351N  1 

SMCAR-CC  1 

SMCAR 'CCP-A  1 

SMCAR-FSA  1 

SMCAR-FSM-E  1 

SMCAR-FSS-D,  BLDG.  94  1 


SMCAR-IMI-I  (STINFO)  BLDG.  59  2 

0 1 CATINNY  ARSENAL.  NJ  07806-5000 

DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ATTN:  SLC8R-DD-T ,  BLDG.  305  1 

ABERDEEN  PROVING  GROUND,  MD  21005-5066 

DIRECTOR 

US  ARMY  MATERIEL  SYSTEMS  ANALYSIS  ACTV 
ATTN:  AMXSY-MP  1 

ABERDEEN  PROVING  GROUND,  MD  21005-5071 

COMMANDER 
HQ,  AMCCOM 

ATTN:  AMSMC- IMP-L  1 

ROCK  ISLAND,  IL  61299-6000 


NO.  OF 
COPIES 

COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SMCRI-ENM  1 

ROCK  ISLAND,  IL  61299-5000 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENGR  ACTV 
ATTN:  AMXIB-P  1 

ROCK  ISLAND,  IL  61299-7260 

COMMANDER 

US  ARMY  TANK-AUTMV  R&D  COMMAND 
ATTN:  AMSTA-DDL  (TECH  LIB)  1 

WARREN.  MI  48397-5000 

COMMANDER 

US  MILITARY  ACADEMY  1 

ATTN:  DEPARTMENT  OF  MECHANICS 
WEST  POINT,  NY  10996-1792 

US  ARMY  MISSILE  COMMAND 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ATTN:  DOCUMENTS  SECT ]  BLDG.  4484 
REDSTONE  ARSENAL,  AL  35898-5241 

COMMANDER 

US  ARMY  FGN  SCIENCE  AND  TECH  CTR 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 

COMMANDER 

US  ARMY  LABCOM 

MATERIALS  TECHNOLOGY  LAB 

ATTN:  SLCMT-IML  (TECH  LIB)  2 

WATERTOWN,  MA  02172-0001 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 

CENTER.  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST  (CONT'D) 


NO.  OF  NO.  OF 

COPIES  COPIES 


COMMANDER 

US  ARMY  LABCOM,  ISA 

ATTN:  SLCIS-IM-TL  1 

2800  POWDER  MILL  ROAD 
ADELPHI,  MD  20783-1145 

COMMANDER 

US  ARMY  RESEARCH  OFFICE 

ATTN:  CHIEF,  IPO  1 

P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709-2211 
DIRECTOR 

US  NAVAL  RESEARCH  LAB 

ATTN:  MATERIALS  SCI  &  TECH  DIVISION  1 

CODE  26-27  (DOC  LIB)  1 

WASHINGTON,  D.C.  20375 


COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MN  1 

EGLIN  AFB ,  FL  32542-5434 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MNF 

EGLIN  AFB,  FL  32542-5434  1 

MIAC/CINDAS 
PURDUE  UNIVERSITY 
2595  YEAGER  ROAD 

WEST  LAFAYETTE,  IN  47905  1 


DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ATTN.  SLCBR-IB-M  (DR.  BRUCE  BURNS)  1 
ABERDEEN  PROVING  GROUND.  MD  21005-5066 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 

CENTER.  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


